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In this talk I summarize recent work done in collaboration with Cliff Burgess and 
Denis Michaud in which we performed a detailed investigation of how solar 
neutrinos propagate through helioseismic waves. We find that the MSW solar 
neutrino spectrum is not modified at all in the presence of seismic waves. This 
finding differs from earlier estimates mainly because most helioseismic waves are 
too weak in the vicinity of the MSW resonance to be of relevance for neutrino 
propagation. A special class of waves may however by subject to an instability and 
potentially have very large amplitudes. These waves do have long wavelengths, a 
situation for which the formalism employed in earlier analyses does not apply. Our 
numerical simulation significantly reduces their influence on neutrino propagation. 



1 Introduction 

The MSwi effective hamiltonian description of neutrino resonant conversion 
in matter provides, by altering the neutrino spectrum, an excellent solution t 
the solar neutrino problem and incidentally the best fit to solar neutrino data 
to date. The promise that the solar neutrino spectrum will be more precisely 
determined in the future spurred investigations of how the MSW mean-field 
description might by altered in the presence of density fluctuations. 

A very efficient mechanism for modifying the MSW spectrum has since 
been proposed uu, It involves the diffuse damping of neutrinos from the co- 
herent MSW evolution as they interact with density fluctuations, and it could 
give rise to significant effects, provided perturbations with correlation length 
very short compared to the neutrino oscillation length and the density scale 
height are actually present in the vicinity of the MSW resonance point. 

What has been missing so far is a plausible source for such delta-correlated 
fluctuations within the sun. We argue that they cannot come about as a su- 
perposition of helioseismic waves, the only plausible source for density fluc- 
tuations deep in the solar interior. Our investigation shows that only a few 
seismic modes can possibly reach amplitudes large enough to be of relevance. 
A random superposition of these modes, however, does give rise to fluctuations 
with, long spatial correlation length — a situation for which the formalism used 
in all breaks down. 
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Here we explicitly compute the wave profiles generated by these seismic 
modes and numerically evolve neutrinos through it, without making use of 
the short-correlation-length approximation. We find that the standard MSW 
neutrino spectrum remains unaltered. 

Below we briefly summarize the approach used in earlier analysesErO. Sec- 
tion 3 then explains why this approach is only valid in the short correlation 
length regime. Section 4 discusses density fluctuations in the sun. Our results 
and conclusions are summarized in section 5. 



2 The Story So Far 

Neutrino flavor evolution in a unperturbed matter background is described, in 
the interaction picture, by 

^ = -i[V Y AC + V M sw(t),p f }, (1) 

where pf is the neutrino flavor density matrix, Vyac ~ k + m ^ L + . . . is the 
vacuum mixing with k and m being the neutrino momentum and mass matrix 
respectively and Vmsw = V / 2Gi?.g e n e (i) denotes the effective interaction with 
the matter background with n e (t) being the electron density, gr e =diag(l,0) the 
charged current coupling matrix and Gf Fermis constant. When integrated 
Eq. (|l|) gives rise to the familiar Parke formulaQ for the MSW survival proba- 
bility. 

In the presence of density fluctuations Eq. ([[]) is modified to becomes 

= -i[VvAC + (Vusw(t)),p f ] (2) 
2G%A(t) [(g e ) 2 Pf + Pf {g e f - 2g e Pf g e ] + 0{G%), 

where (. . .) denotes the average over an ensemble of density fluctuations, so 
that (n e (t)) is the mean electron density. The fluctuations are then contained 
within 5n e = n e — (n e ). The correlation integral A(t) is given by 

A{t) = I (Sn e (t)Sn e (T))dT, (3) 



where t' denotes the initial time of neutrino evolution. When integrated Eq. (||) 
gives rise to the generalized Parke formula for the neutrino survival probability 

p e (m 

P e {t) = \+{\-Pj] Acos20 m (i')cos20 m (*). (4) 
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Here 8 m is the neutrino mixing angle in matter and Pj denotes the standard 
(nonadiabatic) MSW jump probability!!! The damping factor A is given by 



A = exp 



-2G% J A{t) sin 2 20 m {T)dT 



(5) 



Notice that in the absence of fluctuations, Sn e = 0, A(t) vanishes and Eqs. (f|) 
and (0) reduce to the standard Parke formulaB and to Eq. (Q) respectively. 

3 Dealing With Long Correlation Lengths 

The formalism described in the previous section is only valid as long as the 
correlation lengths of the involved density fluctuations is much smaller than 
the neutrino oscillation length. To see this it is useful to expand the electron 
density into a complete set of, respectively uncorrelated, orthonormal func- 
tions: 



n e (t) = (n e (t)} 



(6) 



where the coefficients C n are assumed to be gaussian distributed with vanishing 
mean, (C n ) = 0, and uncorrelated with each other, (C n C m ) — V n 5 nm . These 
basis functions could for example represent helioseismic waves. 

For further illustration it is instructive to consider a one dimensional toy 
model, the 'cell' model, in which the density fluctuation modes <fi n (t) are 
cells of fixed length '£', aligned one after the other so that (6n e (t)6n e (T)) = 
e 2 (n e (i))(n e (r)) if t and r are both within the same cell and zero otherwise. 

For this model the correlation coefficient approximates to A ~ £e 2 n 2 and 
is obviously growing with £. Since the second order term of Eq. (|^) grows 
with A the expansion on which this equation is based will inevitably break 
down at one point. This becomes more transparent when one realizes that the 
involved expansion parameter is not the weak interaction, Gf , but rather a 
combination of weak interaction (Gf), amplitude of the perturbation (e) and 
correlation length (I). 

To get a feeling for when this happens we performed a numerical average 
over the ensemble of density profiles. Since any given neutrino sees only one 
specific density profile as it passes through the sun, we can compute its survival 
probability using the standard MSW evolution (without fluctuations), Eq. (Q), 
taking n e = (n e ) + 5n e . In doing so for a statistical set of density profiles Sn e 
we get an ensemble average for the neutrino survival probability which is valid 
even for long correlation length 0. 

a To save computer time we used the standard Parke formula (without fluctuations) instead 
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Figure 1: Comparing the two averaging procedures in the case of the 'cell' model. The 
neutrino parameters are: ^ 10- 6 Ev 2 /MeV, sin 2 26 v = 0.01 and t 2 = .Q-01. The thin 

solid curve shows the result obtained using the generalized Parke formula Eq. (W) , whereas the 
thick solid curves represents the result of a direct numerical ensemble average. The horizontal 
line represents the standard MSW survival probability, in the absence of fluctuations. 



Where and how the pertrubative formalism, Eqs. (0,^), breaks down is il- 
lustrated in Fig. (0), where we compared the two methods in the framework of 
the 'cell' model described above. The thin solid line shows the result obtained 
with the generalized Parke formula, Eq. (0), whereas the thick solid line rep- 
resents the method described above using a statistical sample of 200 density 
profiles. 

As can be seen from this Figure the numerical ensemble average begins to 
deviate significantly from the generalized Parke formula for I > a few 1000 km. 
For even larger correlation length the only effect of the perturbation is a shift 
of the MSW resonance point, which doesn't affect the survival probability very 
much. Consequently the result of the numerical ensemble average approaches 
the standard MSW result for large correlation length. 

From Fig. (g) we learn two important points: 

• The effect of fluctuations is largest when their correlation length is about 
the length scale at which the perturbative approach breaks down, typically of 

of numerically integrating the full evolution equation Eq. (jlj). The difference turns out to 
be negligible. 
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the order of the neutrino oscillation length. 

• Using the generalized Parke formula beyond its domain of validity leads 
to a severe overestimation of the influence of perturbations on neutrino prop- 
agation. 

4 Helioseismology And Density Fluctuations Inside The Sun 

Our question of concern is what kind of density perturbations in the sun (for 
a good review of the solar interior seeQ) can affect neutrino propagation. An 
inspection of Eq. (^) shows that only density perturbations in the vicinity 
of the MSW resonance point (where sin 2 29 m is maximised) can significantly 
modify the survival probability []. 

Since the MSW resonance typically occurs within the inner 0.5 solar radii, 
density fluctuations in the outer shell of the sun, i.e. the convective zone, and 
the surface won't be of.. relevance. Furthermore a mixed, convective core has 
recently been excluded a. Barring the exotic option of strange hydrodynamics, 
which is not described as perturbations on a spherically symmetric background 
and which must not mix the core, deep inside the sun, this leaves us with 
helioseismic waves as a candidate for perturbations close to the MSW resonance 

P° int - nn 

There are two basic types of helioseismic waves LJ'lI, called 'pressure' (p-) 

waves and 'buoyancy' (g-) waves respectively, which owe their names to the 

main restoring force that acts on a displaced element of matter, being pressure 

and gravity respectively. Other features which distinguish the two types of 

waves are their periods (above/below 30 minutes for g-/p- waves) and the 

regions which they populate. Buoyancy waves typically are strongest in the 

central region of the sun whereas pressure waves are most prominent in the 

convective zone. 

Only pressure modes of low angular degree_penetrate deeply into the solar 
interior. Helioseismic measurements howeverlij indicate that their amplitude 
is far to small to be of relevance to neutrino propagation^. 

This leaves us with buoyancy waves, which typically have their largest 
amplitudes right in the vicinity of the MSW resonance. Since they are damped 
in the gravitationally unstable convective zone they are generically hard to 
detect - and in fact haven't been observed so far - and so could potentially 
be very strong, at least on purely phenomenological grounds. However an 
inspection of the excitation and damping mechanisms indicates that, if at all, 
only a few g-modes with radial wavenumber n < 3 can have large energies, at 

b Even though Eq. ^ only holds for short correlation lengths this remark also proves true 
in the long correlation length regime. 
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Figure 2: Shown is the survival probability for electron neutrinos propagating through a 
& = 13, n = 1 buoyancy wave. The various dashed and dot-dashed lines represent the results 
obtained using the generalized Parke formula with wave amplitudes Sn e /n e as indicated. 
Conversely the solid and dotted line show the result of a direct numerical integration, largely 
indistinguishable from the standard MSW profile. For simplicity of comparison all neutrinos 
have been assumed to be originating from the centre. 

least IO 35 - 37 Ergs This is because these waves could be linearly instable 
and so would grow exponentially until nonlinearities saturate their growthlI3. 
Such modes would correspond to 5n e /n e ~ 10 -7 , potentially giving rise to 
perturbations in the vicinity of 10~ 4 '~ 5 when superimposed — more if they 
were more energetic. 

Since all of these instable modes have long correlation lengths and peak in 
the same area, they will, when superimposed noLeive rise to delta correlated 
white noise as has been assumed in earlier worksffl. The formalism of section 
2 does therefore not apply to them. 

5 Results And Conclusions 

We have computed the solar neutrino spectrum using density profiles gen- 
erated from these linearly instable g-modes. In doing so we averaged over 
neutrino productiop_sites, using the production distributions of the Bahcall- 
Pinsonneault SSMI13, and took into account details such as double resonant 
crossing. 
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We find that the solar neutrino spectrum is not affected on a measurable 
level by helioseismic waves. This is shown in Fig. (Q) where we have evolved 
the neutrinos through a helioseismic g-wave of angular degree I = 13 and radial 
order n = 1. Fig. ^ nicely illustrates the fact that applying the generalized 
Parke formula can lead to a gross overestimate of the size of the effect. The 
direct numerical integration shows no measurable effect at all even for a wave 
amplitude as large as 100%. 

The main reasons of why we find such a small effect are: 

• Most seismic waves are much weaker than necessary to modify neutrino 
propagation. 

• The few waves which could potentially be very strong have very long 
correlation lengths, putting us in the regime where their influence on neutrino 
propagation is significantly reduced (putting us in the right half of Fig. (@)) 
compared to density fluctuations with shorter correlation lengths. 

• Finally, averaging over production cites smears out part of the adia- 
batic dip of the MSW spectrum and provides an important background to any 
fluctuation effect. 

We conclude that density perturbations in the sun are very unlikely to 
modify the MSW neutrino spectrum and the MSW solution to the solar neu- 
trino problem. They certainly do not do so if they are composed of helioseismic 
waves. 
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